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From  a  dataset  of  116  earthquakes,  all  available  SRO  (Special  Research 
Observatory)  and  PSWF  (Preliminary  Signal  Waveform  File)  waveforms  were 
presented  on  an  interactive  computer  graphics  screen  and  then  examined  for 
presence  of  depth  phases  pP,  sP  and  PcP.  Using  the  results,  a  focal  depth  for 
66%  of  the  total  was  postulated.  About  half  of  these,  or  37%,  agreed  with 
depths  published  by  NEP  (Network  Event  Processor)  or  PDE  (Preliminary  Deter¬ 
mination  of  Epicenter)  bulletins.  Earthquakes  of  magnitude  mj,  >  5  are 
likely  to  produce  clear  depth  phases.  Rapid  routine  processing  of  a  large 
number  of  events  by  the  aligned-on-P  method  can  be  achieved  with  interactive 
computer  graphics. 
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INTRODUCTION 


In  the  context  of  monitoring  underground  nuclear  explosions,  the  ability 
to  establish  with  certainty  the  depth  of  an  event  is  an  extremely  powerful 
discriminant.  Although  depth  estimation  can  be  achieved  by  computer  programs 
using  only  first  arrivals,  a  more  accurate  and  reliable  depth  estimation  can 
be  made  if  the  depth  phases  pP  and  sP  are  observed.  Body  waves  reflected 
back  from  the  surface  of  the  earth  comprise  pP  or  sP,  and  similar  reflections 
from  the  core  are  called  PcP.  These  phases  arrive  within  2  minutes  after  P, 
and  are  imbedded  in  the  P-wave  coda. 

The  primary  factors  that  influence  the  appearance  (or  non-appearance)  of 
a  depth  phase  are  thought  to  be: 

magnitude  of  earthquake 
depth 

focal  mechanism 

source-to-receiver  orientation 
signal  to  noise  ratio  at  receiver 
structure  around  the  hypocenter. 

The  seismological  literature,  to  the  author’s  knowledge,  contains  no 
information  on  the  observability  of  depth  phase  and  no  work  that  systemati¬ 
cally  relates  the  above  factors  to  actual  observations.  However,  Dahlman's 
and  Israelson’s  (1977)  book  includes  a  chapter  on  depth  estimation.  They  cit 

two  reports,  one  by  Lacoss  (1969)  and  the  other  by  Yamamoto  (1974),  that  haw 
examined  and  reported  on  the  frequency  of  observation  of  depth  phases.  The 
location  of  events  in  the  Lacoss  study  is  defined  by  the  LASA  beam,  while 


Dahlman,  0.  and  H.  Israelson,  1977,  Monitoring  Underground  Nuclear  Explosions 
New  York,  Elsevier  Scientific  Publishing  Co. 

Lacoss,  R.  T.,  1969,  A  large  population  LASA  discrimination  experiment :  MIT 
Lincoln  Laboratory,  Technical  Note  1969-24. 

Yamamoto,  M. ,  1974,  Estimation  of  focal  depth  by  pP  and  sP  phases:  National 
Defense  Research  Institute,  Stockholm,  Report  A4505-A1. 


Yamamoto's  work  was  concerned  with  Japanese  earthquakes.  Chiburis  and  Ahner 
(1969)  produced  a  third  key  report  covering  earthquakes  in  North  and  South 
America. 

In  our  study,  waveforms  from  a  global  network  of  seismometers  are  pre¬ 
sented  on  a  computer  graphics  screen  and  then  examined  by  an  analyst  for 
evidence  of  depth  phases.  The  goal  is  to  obtain  information  on  how  frequently 
a  sparse  global  network  of  seismic  stations  can  observe  earthquake  depth  phases. 
The  data  used  consists  of  116  events  derived  chiefly  from  the  Network  Event 
Processor  (NEP)  bulletin.  Earthquake  locations  are  worldwide,  rather  than 
limited  to  specific  regions  as  in  earlier  studies. 

While  the  three  earlier  studies,  as  well  as  this  study,  differ  in 
fundamental  design  (e.g.,  Yamamoto's  study  did  not  utilize  waveform  data), 
they  agree  broadly  in  their  conclusions.  That  is,  that  depth  phases  can  be 
observed  in  about  two  thirds  of  all  m^  a  5  events,  and  that  half  the 
implied  depths  agree  with  USGS  depths.  The  major  difficulty  in  observing  depth 
phases  is  identifying  a  phase  imbedded  in  the  P-coda,  particularly  if 
its  amplitude  is  equal  to  or  less  than  the  coda  amplitude  at  the  time  the 
depth  phase  appears.  Presenting  available  waveforms  on  a  graphics  screen, 
together  with  an  interactive-predictor  mark,  helped  the  analyst  make  a 
decision. 

Results  from  this  project  will  be  further  analyzed  with  a  computer  program 
utilizing  Pearce's  (1977)  ideas  for  establishing  fault  plane  solutions  using 
P-pP  and  P-sP  amplitude  ratios.  While  depth  phases  are  not  always  observed  at 
all  network  stations,  their  absence  is  as  important  as  their  presence  because 
they  act  as  a  constraint  on  possible  fault  solutions. 


Chiburis,  E.  F.  and  R.  0.  Ahner,  1969,  The  comparative  detectability  of  pP 
at  LASA,  TSFO,  UBSO  and  CPSO:  Seismic  Data  Laboratory  Report  No.  231, 
Teledyne  Geotech,  Alexandria,  Virginia. 

Pearce,  R.  G.,  1977,  Fault  plane  solutions  using  relative  amplitudes  of  P  and 
pP:  Geophys,  J.  R.  Astr.  Soc.,  50,  381-394. 
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DETECTION  OF  DEPTH  USING  COMPUTER  GRAPHICS 


One  method  of  enhancing  the  capability  to  detect  and  identify  depth  phases 
pP,  sP  and  Pc P  is  to  arrange  the  short  period  signals  from  a  network  of 
seismometers  so  they  are  aligned  on  first  motion  of  P  onset.  Presumably, 
secondary  depth  arrivals,  if  present,  will  form  a  coherent  pattern  with 
evidence  of  moveout  of  the  depth  phase  with  increasing  epicentral  distance, 
although  the  moveout  may  be  slight.  The  moveout  is  slight  because  the  time 
delays  between  the  arrival  of  pP  and  P,  as  well  as  sT  and  P,  are  chiefly 
functions  of  focal  depth.  They  are  almost  independent  of  epicentral  distances 
for  events  less  than  100  km  deep,  so  the  moveout  is  small,  perhaps  just  2  or 
3  seconds,  according  to  the  range  of  epicentral  distances  covered  by  the 
available  signals.  However,  producing  this  type  of  section  with  hard  copy  is 
time-consuming  and  not  suitable  for  routine  event  analysis.  This  report 
describes  an  investigation  into  using  computer  graphics  to  build  up  the 
aligned-on-P  section  for  116  events,  and  expresses  the  results  in  terms  of 
how  frequently  depth  phases  are  seen,  how  they  are  identified  and  the  level 
of  confidence,  as  well  as  the  implications  for  on-line  discrimination  between 
explosions  and  earthquakes. 

Ihe  experiment  was  designed  to  use  the  PDP-15  computer  to  acquire  all 
signals  available  at  the  SDAC  for  events  listed  in  the  NEP  bulletin  for 
the  period  1  December  1977  to  15  December  1977,  to  then  present  these  signals 
on  the  graphics  screen  so  they  could  be  aligned  on  P-onset ,  and  to  give  the 
analyst  a  moveable  vector  for  pP  calculated  arrivals.  The  signals  are  pre¬ 
sented  in  order  of  increasing  epicentral  distance.  The  analyst  shifts  each 
signal  left  or  right  to  align  P-onsets  to  a  vertical  reference  line,  then 
enters  a  depth  from  teletype  keyboard,  and  a  series  of  vectors  appear  on  the 
screen  connecting  the  baselines  of  each  signal  at  the  point  corresponding  to 
the  calculated  arrival  position  of  the  pP  phase.  Also  shown  on  the  screen 
are  the  arrival  positions  for  sP  and  PcP  for  the  chosen  depth.  Computer 
response  is  essentially  instantaneous  for  entry  and  display  of  a  new  depth 
vector. 

Two  data  bases  were  used  to  acquire  and  present  signals:  the  SRO  network 
and  the  PSWF  tape  from  NEP.  At  the  time  of  the  survey  the  SRO  network  consisted 


of  11  reporting  stations.  Three  Alaskan  stations,  plus  LASA  from  the 
PSWF  tape,  gave  a  total  of  15  potential  signals  for  the  display.  However, 
the  maximum  number  of  signals  seen  for  any  one  event  was  9.  The  PDP-15 
graphics  screen  can  hold  11  signals.  The  SDAC  bulletin  reported  a  total  of 
116  events  in  the  period  1  through  15  December  1977.  Using  the  computer 
graphics  scheme,  we  observed  depth  phases  of  sufficient  clarity  to  postulate 
a  focal  depth  for  77  of  these  events,  or  66%  of  the  total  number  of  bulletin 
events. 

The  116  events  produced  a  total  number  of  658  signals  at  the  network's 
reporting  stations.  Of  these  signals  the  analyst  picked  312  signals,  or  ■  50% 
as  containing  a  depth  phase.  The  pick  was  not  made  on  an  individual  signal, 
but  rather  it  was  made  while  the  graphics  screen  presented  all  signals  de¬ 
tected  for  that  event,  even  though  some  of  the  signals  there  did  not  contain 
a  depth  phase,  or  might  contain  conflicting  and  extraneous  phases. 

The  analyst  recorded  a  confidence  rating  to  the  depth  assignment  ranging 
from  1  =  (high  confidence,  numerous  clear  depth  phases)  to  4  =  (unable  to 
assign  a  focal  depth).  This  rating  was  subjective;  it  depended  upon  the  number 
of  stations  reporting,  fraction  of  signals  on  the  screen  showing  a  depth  phase, 
and  non-ambiguity  of  the  ptiase.  The  number  of  events  and  percentages  of  the 
total  number  in  each  of  these  categories  are: 


Confidence  in 

Depth  Assignment  (1  =  High) 

Number  of 

Events 

Percentage  of  Total 

1 

12 

10% 

2 

31 

27% 

3 

34 

29% 

4  (no  depth  assign¬ 
ment  possible) 

39 

34% 

116 

100% 

The  66%  rate  of  depth  postulation  compared  reasonably  with  the  study  of 
150  Eurasia  events  that  Lacoss  (1969)  recorded  at  LASA.  Lacoss  found  that  a 
secondary  phase  could  be  seen  in  about  60%  of  the  events,  and  that  of  these. 
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about  80%  could  be  associated  with  pP  or  sP.  Thus,  correctly  identifying  a 
deptli  phase  was  possible  in  about  50%  of  the  events  studied  at  LASA. 


However,  where  only  a  single  receiving  station  exists,  such  as  the  LASA 
beam,  establ isliing  whether  a  solitary  secondary  phase  is  pP  or  sP  is  virtually 
impossible  if  other  information  is  unavailable.  Still,  it  is  statis¬ 
tically  sound  to  consider  an  unassigned  depth  phase  to  be  pP,  because 
Yamamoto  (1969)  showed  that  for  events  in  the  Japan  region  the  occurrence  of 
pP  was  six  times  as  great  as  the  occurrence  of  sP. 

Table  1  shows  that  this  method's  applicability  and  the  confidence  level 
of  a  depth  assignment  are  strong  functions  of  the  magnitude  of  the  116  events. 
Only  40%  of  events  with  m^  less  than  5.0  could  be  assigned  a  focal  depth, 
while  83%  of  events  with  m^  greater  than  or  equal  to  5.0  showed  sufficient 
evidence  of  depth  phases  to  be  assigned  a  focal  depth.  In  addition,  events 
with  smaller  magnitudes  tended  to  cluster  at  the  lower  confidence  levels  of 
depth  assignment,  while  the  greater  magnitude  events  tended  toward  a  high 
confidence  level  of  1  in  their  depth  assignment. 

Depth  Verification 

Verification  of  the  focal  depths  assigned  by  the  analyst  based  on  his 
interpretation  of  the  information  presented  on  the  computer  graphics  screen 
is  complicated  by  a  number  of  factors.  The  NEIS  depths  do  not  necessarily 
agree  with  the  NEP  depths.  Both  NEP  and  NEIS  compute  focal  depths  in  two 
ways.  Restrained  depth  solutions  (R)  indicate  that  the  NEP  analyst  fixed  the 
depth  during  the  reiterations  of  the  location  program  based  on  additional 
information  such  as  depth  phases.  The  (R)  solution,  designated  (D)  in  the 
NEIS  bulletin,  is  considered  to  be  more  accurate  than  a  free-running  HYPO 
location  solution,  which  is  marked  (F)  in  the  NEP  bulletin,  and  is  without 
added  notation  in  the  NEIS  bulletin. 

In  some  cases  the  graphics  screen  showed  that  more  than  one  focal  depth 
interpretation  could  be  made  for  a  set  of  seismic  traces.  In  these  cases  we 
attempted  to  use  the  depth  with  the  greatest  confidence,  usually  containing 
the  most  traces  which  showed  evidence  of  a  depth  phase. 


TABLE 


Finally,  the  NEIS  bulletin  contains  a  third  category  of  reported 
depths  labeled  33N,  which  indicates  that  it  was  held  at  33  kjn  (normal  depth), 
if  at  any  point  in  the  computation  the  depth  became  negative,  or  was  other¬ 
wise  unsatisfactory,  and  the  earthquake  probably  had  a  shallow  focus. 

In  Figure  1,  the  computer  graphics  depths  are  plotted  against  NEP  depths 
for  those  events  in  the  NEP  bulletin  with  non-zero  focal  depth.  These  NE!’ 
depths  contain  both  free  and  restrained  depths,  shown  by  different  symbols 
on  tile  graph.  The  agreement  with  the  restrained  NEP  depths  is  excellent  and 
is  considerably  better  than  the  match  witli  the  unrestrained  location  solutions. 
All  of  the  computer  graphics  depths  which  differ  witii  NEP  by  more  than  30% 
have  been  interpreted  as  shallow  focus,  70  km  or  less.  None  of  non-agreeing 
events  have  been  interpreted  as  deep  focus  events. 

Figure  2  shows  a  presentation  similar  to  Figure  1  for  NEIS  events  which 
also  have  a  non-33K  depth.  There  is  only  one  restrained  depth  in  this  set, 
agreeing  well  witii  the  computer  graphics  depth.  This  NEIS  data  set  seems  to 
show  more  scatter  than  the  NEP  data  set;  the  explanation  for  this  is  not 
obvious . 

Figure  3  is  a  histogram  oi  computer  graphics  depths  for  28  earthquakes 
reported  as  "33N"  in  the  NEIS  bulletin.  It  is  significant  that  all  the 
computer  graphics  depths  are  shallow  (<_  7U  km),  and  that  the  peak  oi  the 
distribution  occurs  close  to  20  km,  thus  tending  to  support  the  NEIS  procedure 
of  assigning  3  IN  when  the  location  program  does  produce  a  satisfactory  depth, 
and  it  is  a  geophysicist's  opinion  that  the  focus  is  probably  shallow. 

Nineteen  of  these  twenty-eight  earthquakes  are  contained  in  the  NEP  bulletin 
where  they  are  restrained  to  surface  focus,  indicating  also  that  the  NEP 
computer  location  program  produced  unsatisfactory  focal  depths.  There! ore, 
it  is  possible  that  the  computer  graphics  approach  has  succeeded  in  assigning 
focal  depths  to  shallow  earthquakes  when  location  programs  using  arrival 
times  fail. 

Unfortunately,  there  is  no  method  readily  available  to  verify  these 
depth  determinat  ions  in  the  cases  where  a  location  program  fails  to  achieve 
a  satisfactory  depth.  It  is  not  surprising  that  stacking  a  number  of 
aligned-on-P  waveforms  should  help  in  recognizing  any  coherent  pattern  that 
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Figure  1.  Depth  determined  with  the  aid  of  a  computer 
graphics  screen  presentation  of  all  P  phases 
plotted  against  the  depth  for  the  same  event 
in  the  N'EP  bulletin. 


may  be  presented,  and  the  method  would  probably  be  more  effective  If  more 
signals  were  available. 

To  get  an  overall  picture  of  the  performance  of  this  project,  we  now 
look  at  how  many  depth  determinations  agree  or  disagree  with  NEP  and  NEIS. 
Lines  corresponding  to  plus  and  minus  50%  of  bulletin  depth  are  shown  on 
Figures  1  and  2;  these  lines  are  used  to  arbitrarily  class  the  events  into 
those  that  are  "verified"  and  those  that  are  not.  Furthermore,  we  make  the 
assumption  that  all  the  "33N"  events  are  "verified." 

Counting  in  this  way  yields  the  following  results: 


Agree  with  NEP  28 
Agree  with  NEIS  25 
Agree  with  "33N"  28 
No  agreement  27 

Total  102 


Of  course,  this  totals  more  events  than  the  77  events  that  were  assigned 
a  focal  depth,  because  if  an  event  is  on  both  NEIS  and  NEP  bulletins,  it  is 
counted  twice  in  the  above  tabulation.  However,  the  percentage,  here  80%,  of 
computer  assisted  focal  depth  determinations  that  have  been  "verified"  is 
correct. 
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DISCUSSION  AND  CONCLUSIONS 

Table  II  compares  the  results  of  this  study  with  the  studies  of  Lacoss, 
Yamamoto,  and  Chiburis  and  Ahner  from  the  point  of  view  of  the  number  of 
earthquakes  that  have  observable  depth  phases. 

Chiburis  and  Ahner  (1969)  suggest  that  the  large  array  is  a  better 
detector  of  pP  and  sP,  but  the  information  in  a  single  beam  is  insufficient 
to  choose  between  these  two  phases.  They  reached  this  conclusion  in  their 
study  of  the  comparative  detectability  of  later  seismic  phase  arrivals  of 
LASA  and  a  small  continental  network  of  three  observatories.  However, 
several  factors  appearing  in  Table  II  make  a  comparative  interpretation  of 
these  four  studies  difficult.  These  are:  (1)  a  different  geographical 
region  is  sampled  in  each  case;  (2)  Chiburis  and  Ahner' s  data  base  must 
have  been  selected  to  eliminate  events  that  did  not  show  secondary  arrivals, 
but  they  did  not  report  how  many  events  were  discarded;  and  (3)  Yamamoto’s 
relatively  low  (32%)  fraction  of  total  number  showing  secondary  arrivals 
stems  from  his  requirement  that  at  least  5  stations  (out  of  109)  reported 
later  phases. 

An  important  seismological  question  raised  in  this  report  focuses  on  what 
percentage  of  earthquakes  might  show  later  depth  phases.  A  later  phase  is 
easier  to  recognize  in  deep  (>  100  km)  earthquakes;  the  phases  are  better  sepa¬ 
rated  because  the  coda  are  generally  more  simple  and  die  out  more  quickly  than 
shallow  earthquake  coda.  Both  this  report  and  Yamamoto's  study  show  that  for 
earthquakes  with  m^  >  5.0  the  frequency  of  appearance  of  later  depth  phases  is 
quite  large;  Yamamoto  asserts  that  his  method  has  an  applicability  of  70%,  if 
is  greater  than  5.0.  This  study  found  evidence  of  depth  phases  in  83%  of 
events  where  m^  was  greater  than  5.0.  Thus,  3  out  of  4  earthquakes  of  magnitude 
5.0  or  larger  may  be  expected  to  show  a  later  depth  phase.  Presumably,  the 
same  ratio  persists  for  earthquakes  of  smaller  magnitude,  but  the  number  actu¬ 
ally  observed  is  less  because  of  few  reporting  stations,  less  azimuthal  coverage, 
and  lower  signal  to  noise  ratio.  An  array,  by  itself,  cannot  verify  if  a  phase 
is  pP  or  sP,  but  it  is  more  likely  that  a  phase  is  pP  rather  than  sP  by  about 
6  to  1  (Yamamoto,  1969). 
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Examples  of  Aligned  Seismograms 

The  method  employed  in  this  study  is  as  follows.  Two  tapes  are  mounted 
on  the  PDP-15  tape  drives;  the  PSWF  tape  contains  the  NEP  signals,  and  the  SRO 
Day  tape  contains  all  reporting  stations.  For  individual  events,  arrival  times 
are  computed  and  the  SRO  tape  is  searched;  any  signals  present  are  windowed 
and  merged  with  signals  on  the  PSWF.  The  analyst  studies  the  display  and  makes 
the  best  choice  of  a  focal  depth  based  upon  the  displayed  data. 

Some  examples  of  the  hard  copy  of  the  graphics  screen  are  presented  in 
the  Figures  4  through  17.  The  upper  number  at  the  start  of  a  seismogram  is 
the  epicentral  distance  in  degrees,  the  lower  smaller  number  is  the  azimuth 
in  degrees.  Station  names  are  not  shown  on  the  graphics  screen,  but  are 
available  to  the  analyst  on  the  print  sheet  accompanying  each  event. 
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Depth  =  1*2  km.  42  (NEP) ,  33  (PDE).  Good  depth 


3481252 


3481252,  14  Dec.  I1 
Good  phase  clarity 


3472130 


Event  3395929,  5  Dec.  1970,  Tonga  Islands  Reg 


Figure  12.  Event  3396551,  5  Dec.  1970,  Near  coast  of  Central  Chile 
mj,  =  5.2.  Depth  =  30  km,  0  (NEP) ,  32D  (PDE) .  Good 
agreement  with  PDE  depth  although  phase  quality  is  only 
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3471297 


Event  3471297,  13  Dec.  1977,  Fiji  Islands  Region,  mb 
5.3.  Depth  =  534  (PDE) .  The  display  screen  holds  50 
seconds  of  signal.  If  the  focal  depth  is  deeper  than 
200  km  a  depth  phase  is  off  the  screen.  Deep  events 
show  simple,  explosion  like  P  waveforms. 
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APPENDIX 


Time  delays  of  surface  reflections  (pP  and  sP)  and  the  core  reflection 
PcP,  relative  to  the  direct  P  wave  as  functions  of  epicentral  distance.  A, 
and  focal  depth,  according  to  the  travel-time  table  of  Jeffrey  &  Bullen. 

The  band  for  PcP  covers  the  depth  interval  of  0  -  100  km. 

DEPTH  ESTIMATION 

TRAVEL  TIME 
DIFFERENCE 
(MINUTES) 


Figure  A-l.  Time  delays  of  surface  reflections  (pP  and  sP)  and  the 
core  reflection,  PcP,  relative  to  the  direct  P  wave  as 
functions  of  epicentral  distance.  A,  and  focal  depth, 
according  to  the  travel— time  table  of  Jeffreys  &  Bullen 
(1967).  The  band  for  PcP  covers  the  depth  interval  of 
0-100  km.  (After  Yamamoto,  1974.) 
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